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Abstract 

An evaluation of slow neutron deceleration through their interaction with 
nanoclusters in liquid helium is performed. It is shown that this process is 
strongly suppressed if clusters are bound by the van der Waals interaction. 
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Ultracold neutrons (UCN) HI is now useful instrument in fundamental neutron physics BUO. 
UCN have energies below ~ 300 neV, velocities below ~ 7 — 8 m/s, their characteristic temperature 
is in the mK range. They can be stored in material and magnetic traps and contained during long time 
- very close to the neutron decay lifetime. 

Recently a method of production of high density UCN gas has been proposed in Q. It consists 
in deep cooling of a beam of very cold neutrons (velocities 50-100 m/s) through their collisions 
with nanoparticles (clusters) made from materials with low neutron absorption (deuterium, heavy 
water ice, oxygen, etc), immersed into very cold superfluid helium. The scattering cross section of 
the incident large wavelength neutrons with particles of 1 — 10 nm size is rather large, and if the 
nanoparticles are free, the energy of incident neutrons is reduced due to the energy transfer to the 
recoiling nanoparticles. The calculations of the very cold neutron deceleration and thermalization 
in such a process of scattering by free nanoparticles are described in 0. The test experiments on 
scattering of slow neutrons by deuterium and heavy water ice clusters in liquid helium at a temperature 
1 .6 K have been performed in [J6J. 

1. But it is known that due to the van der Waals interaction the particles are attracted. A sim- 
ple estimate based on the force of the van der Waals attraction [10J gives that at a distance between 
nanoparticles of ~ 100 nm the time needed for particles to come in contact is of the order of mi- 
croseconds. Thus, in the very short time the nanoclusters agglomerate inside the superfluid helium 
into larger clusters and finally form loosely bound highly porous medium - gel consisting of clumps 
connected by strands Gl[8l|9l[6]|. 

Therefore from the very beginning we should consider the slow neutron interaction with bound 
clusters. In the limit of very low temperature of the scattering medium compared to the incident 
neutron energy E n : E n ^> kT we can use the model of neutron scattering by the massive harmonic 
oscillator at T = 0. 

The fraction of the elastic neutron scattering (the ratio of the elastic scattering to the total one) by 
the oscillator in the lattice of similar clusters is determined by the Debye-Waller factor: 

B D _ W = exp(-q 2 < u 2 >), (1) 
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where q is the neutron wave vector transfer, q = 4n sin(9 / '2) / 'A, A is the neutron wavelength, 9 is 
the scattering angle, and < u 2 > is the time averaged mean squared displacement of the oscillator due 
to thermal motion. Generally < u 2 >=< E > jMuj 2 , where < E > is the mean excitation energy of 
the oscillator, M is the mass of the oscillator, uj is the oscillator frequency. At T=0 the displacement 
is determined by the zero oscillations: < u 2 >= %/2Mu. 

Take typical figures for the proposed method [ffl[51|6l: the incident neutron velocity 50m/s (the 
energy ~ 13 aeV , the wavelength A = 8 nm), the scattering angle 9 = ir/3, the scattering wave vector 
q = 0.8 nm -1 , the cluster radius 2 nm, the cluster mass (deuterium) is M ~ 6 X 1CT 21 g ~ 4 x 10 3 
neutron masses. 

At the exponent in the Debye-Waller factor about 1 (the fraction of elastic scattering is about 0.3) 
we need < u 2 > 1 ' 2 = l/q = 1.2 nm, the oscillation energy huj ~ 3.5 neV. 

What is the mean squared displacement and the oscillation energy of these clusters bound in gel? 
Cluster-cluster interaction involves mainly the pair van der Waals forces between molecules belonging 
to different clusters. This pair intermolecular interaction can be described by the Lennard- Jones 
potential: 



U L .j(r) = 4e 



CT\ 12 f O 
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where e is the value of the potential energy minimum, and 2 1 / 6 <r is the position of this minimum. 

For the approximate estimate we can use an isotropic part of the intermolecular deuterium interac- 
tion: e pa 3.2 meV, a w 0.3 nm [fTTT|. 

Cluster-cluster interaction can be calculated by the integration of the pair intermolecular interac- 
tions over the volumes of the clusters. The expressions for these integrals in the case of spherical 
clusters have been obtained in [|T2l[T3l . Using these expressions we can write the expression for the 
interaction of two clusters with the radius R 



n 2 ea s R 

U(r) 



3vi(r-2R) 



a 6 



- 1 



(3) 



_210(r - 2R) 

where v is the volume of the molecule. This expression is valid for closely placed clusters when 
r ~ 2R. 

Oscillation frequency of the cluster in the minimum of this potential uj = (/x/M) 1 / 2 , where 
H = Jj"(r = r min ). 

Taking t> ~ 3.3 x 10~ 2 nm 3 (the deuterium density p pa 0.2 g /cm 3 ) we obtain uj pa 6.3 x 10 11 s~ l , 
the energy Huj pa 0.4 meV, < u 2 >pa 1.5 x 10~ 5 nm 2 . Finally, for the considered above typical slow 
neutron-deuterium cluster scattering the exponent in the Debye-Waller factor < u 2 > q 2 ~ 1 x 10~ 5 
- the probability of the inelastic neutron-cluster scattering is of the order of 10~ 5 . If the clusters are 
not spherical but faceted the interaction could be larger. 

The neutron deceleration by the excitation of vibrational modes of clusters in a lattice takes place 
in the meV neutron energy range, the total neutron-cluster scattering cross section is decreasing as 
an inverse energy. At the neutron energy ~ 20 meV - the depth of the potential well (Eq. (3)) the 
probability of the inelastic scattering is still about 10 -2 . 

It was found experimentally [|71 [H |9l that the clusters in liquid helium are coated with a thin layer 
of solidified helium. These layers screen the nanoparticles from one another and should reduce the 
interaction between them. 

Contrary to the Lennard- Jones interaction between molecules, the van der Waals interaction be- 
tween clusters is decreasing rather slowly with increasing distance between them: it falls only ~ 
three-four times when the gap between deuterium clusters is about 1 nm. Thus it does not seem that 
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the motion of clusters in the solid helium matrix is free enough for the effective energy transfer from 
very slow neutrons to clusters. 

It is not known if it is possible to prevent coagulation in liquid helium. 

2. Another possibility mentioned for cooling of the cold neutron beam down to the UCN energy 
range in the neutron interaction with clusters at T = is the excitation of the internal elastic modes 
of the clusters. 

The theoretical consideration of slow neutron scattering by small particles has a rather long history 
(see for example |[T5lO . The calculations of neutron scattering with excitations of oscillations of bub- 
bles in liquid helium have been performed recently by V. Gudkov [TToTl . His results may be applied to 
slow neutron scattering by clusters if to replace the characteristic frequencies of bubble excitations by 
the frequencies of the elastic modes of clusters. For our example of the deuterium clusters at T = we 
should consider the low lying compressional spheroidal modes with I = (radial breathing modes), 
the spheroidal I = 2 (surface) vibrations and the torsional / = 2 modes. Their eigenfrequencies may 
be calculated from the boundary condition equation for the / = radial modes ifTTll . for the spheroidal 
(surface) vibrations and torsional modes we can use the results from the paper lfT8l on low-frequency 
elastic modes of spherical particle. Using the data for the sound velocities in solid deuterium from 
iTmi : the longitudinal velocity V{ ~ 1.7 x 10 5 cm/ ' s, the transverse velocity v t ~ 1. x 10 5 cm/ ' s, and 
the solid deuterium density p ps 0.2 g/cm 3 we obtain the energy of the first radial / = compressional 
oscillation level of m 1.3 meV. The oscillation energy of the I = 2 spheroidal and torsional modes 
HI is « 1 meV. 

We rewrite the Eq. (13) of the publication [fT6ll for the neutron down-scattering cross sections with 
excitation of the first internal modes in somewhat different form: 



where A is the ratio of the cluster mass to the neutron mass, V is the Fermi potential of the cluster, 
E n and E int are the energies of the neutron and of the internal mode of the cluster, respectively, R 
is the radius of the cluster, £ = qR, £ min and ^ max are the minimum and maximum ^-values in the 
neutron scattering by the cluster with the energy transfer from the neutron to the cluster AE n = E int , 
and $(£) is the corresponding mode form factor |[T6ll . 

Taking A = 4 x 10 3 , R = 2 nm, V = 100 neV, E n = 2 meV and E int = 1 meV we obtain that 
the corresponding cross sections are very small: for the excitation of I = radial mode 07=0 ~ 17 mb, 
for the excitation of the spheroidal surface I = 2 mode it is even lower: cr i=2 ~ 15 fib. 
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